INTERACTIONS OF MASSIVE STARS WITH THEIR PARENTAL CLOUDS 



On 
On 



C/3 



> 
O 



Jose Franco 1 , Guillermo Garcfa-Segura 1 & Tomasz Plewa 2 
ABSTRACT 

Here we discuss the interaction of massive stars with their parental molecular clouds. 
A summary of the dynamical evolution of HII regions and wind-driven bubbles in high- 
pressure cloud cores is given. Both ultracompact HII regions and ultracompact wind- 
driven bubbles can reach pressure equilibrium with their surrounding medium. The 
structures stall their expansion and become static and, as long as the ionization sources 
and the ambient densities remain about constant, the resulting regions are stable and 
Oh 1 long lived. For cases with negative density gradients, and depending on the density 

distribution, some regions never reach the static equilibrium condition. For power- 
law density stratifications, p oc r~ w , the properties of the evolution depend on a 
critical exponent, w cr it, above which the ionization front cannot be slowed down by 
\ recombinations or new ionizations, and the cloud becomes fully ionized. This critical 

exponent is w cr n = 3/2 during the expansion phase. For w > 3/2 the gas expands 
supersonically into the surrounding ionized medium, and there are two regimes separated 
by w — 3. For 3/2 < w < 3, the slow regime, the inner region drives a weak shock moving 
with almost constant velocity through the cloud. For w > 3, the fast regime, the shock 
becomes strong and accelerates with time. Finally, the evolution of slow winds in highly 
. pressurized region is described briefly. 

O 

Q\. 1. INTRODUCTION 

^3 

Young stars display vigorous activity and their energy output stirs and heats the gas in their vicinity. Low- 
mass stars provide a small energy rate and affect only small volumes, but their collective action can provide 
partial support against the collapse of their parental clouds, and could regulate some aspects of the cloud 
evolution (e.g., Norman & Silk 1980; Franco & Cox 1983; Franco 1984; McKee 1989). In contrast, stars with 
initial masses above 8 Mq, massive stars, inject large amounts of radiative and mechanical energy from their 
moment of birth until their final explosion as a supernova. In the general, low-density, interstellar medium of a 
gaseous galaxy, the combined effects of supernovae, stellar winds, and HII region expansion destroy star-forming 
clouds, produce the hottest gas phases, create large expanding bubbles, and are probably responsible for both 
stimulating and shutting off the star formation process at different scales (e.g., Cox & Smith 1974; Salpeter 
1976; McKee & Ostriker 1977; Franco & Shore 1984; Cioffi & Shull 1991; Franco et al. 1994; Silich et al. 
1996). Thus, the collection of OB associations in gaseous galaxies represents a rich energy source which may be 
controlling the general structure of the interstellar medium, and the star formation rate (e.g., Mueller & Arnett 
1976; Gerola & Seiden 1978; Franco & Shore 1984; Dopita 1988; see reviews by Tenorio-Tagle & Bodenheimer 
1988, Franco 1991,1992, Ferrini 1992, and Shore k Ferrini 1994). 

The strong UV radiation field from massive stars create large photoionized, HII regions. Young HII regions 
have large pressures and, when the pressure of the surrounding medium is low, they expand fast and drive a 
strong shock wave ahead of the ionization front. Expanding HII regions ionize and stir the parental cloud and, 
when the ionization front encounters a strong negative density gradient, they create fast "champagne" flows 
and can also generate comctary globules and elephant trunks (e.g., Tenorio-Tagle 1982; Yorke 1986; Franco et 
al. 1989,1990; Rodriguez-Caspar et al. 1995; Garcia-Segura & Franco 1996). These flows are generated by the 
pressure difference between the HII region and the ambient medium, and they are responsible for the disruption 
of the cloud environment (e.g., Whitworth 1979; Elmegreen 1983; Larson 1992; Franco et al. 1994). Individual 
HII regions are bright objects and they are used as tracers of the active star formation sites in external galaxies 
(e.g., Osterbrock 1989). 
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Stellar winds and supernova explosions, on the other hand, being powerful sources of mechanical energy 
generate overpressured regions which drive shock waves into the ambient medium (see reviews by Ostriker & 
McKee 1988, and Bisnovatyi-Kogan & Silich 1995). The resulting wind-driven bubbles and supernova remnants, 
either from a single progenitor or from an entire association, are believed to generate most of the structuring 
observed in a gaseous galactic disk (e.g., Reynolds & Ogden 1978; Cowie et al. 1979; Heiles 1979; McCray & 
Snow 1979). Actually, many observed structures in the Milky Way and in external galaxies have been ascribed 
to this stellar energy injection (e.g., Heiles 1979, 1984; Brinks & Bajaja 1986; Dcul & Hartog 1990; Palous et 
al. 1990, 1994). These bubbles may create fountains or winds at galactic scales (e.g., Shapiro & Field 1976; 
Chevalier & Oegerle 1979; Bregman 1980; Cox 1981; Heiles 1990; Houck & Bregman 1990), and their expanding 
shocks have also been suspected of inducing star formation (e.g., Hcrbst & Assousa 1977; Dopita et al. 1985). 
Thus, stellar activity creates a collection of cavities with different sizes, and can be viewed as an important 
element in defining the structure and activity of star-forming galaxies. Galaxies, however, are open systems 
and their properties are also defined by the interactions with neighboring galaxies. Here we describe only the 
effects of the stellar energy injection in high-density, high-pressure, regions. 

2. THE PARENTAL CLOUDS AND MASSIVE STARS 

Molecular clouds have complex density and velocity distributions, and are composed of a variety of high- 
density condensations. In our Galaxy, they have non-thermal turbulent velocities, reaching up to about 10 
km s _1 , and fairly strong magnetic fields, of up to tens of mG (e.g., Myers & Goodman 1988). The average 
densities for molecular cloud complexes is between 10 2 and 10 3 cm~ 3 , but the high-density condensations have 
average densities of about ~ 10 6 cm -3 and may even reach values in excess of 10 8 cm -3 (e.g., Bergin et al. 
1996; Akeson et al. 1996; see review by Walmsley 1995). These high density condensations, or cloud cores, 
are the actual sites of star formation, and the initial shape and early evolution of the resulting HII regions 
depend on the corresponding core density distributions and pressures. Theoretical studies on the collapse of 
clouds indicate that rotating and magnetized cores evolve into flattened (disk-like) structures (e.g., Bodenheimer 
& Black 1978; Cassen et al. 1985), but nonrotating and nonmagnetic cases remain spherically symmetric with 
power-law density distributions (e.g., Larson 1974). Isothermal spheres in hydrostatic equilibrium have a density 
distribution p <~ r~ 2 , and the distribution evolves towards r~ 3 / 2 during the free-fall collapse. The pressures at 
the centers of these cores are fairly large, and can reach values of about 5 or 6 orders of magnitude above the 
pressures at the cloud boundaries (see Garcfa-Segura & Franco 1996). 

Recent observational studies are revealing the structure of the dense star forming regions (see review by 
Walmsley 1995), and provide the average parameters of the core density distributions. Radio observations of 
cloud cores and dark clouds, along with visual extinction studies in nearby star forming clouds, indicate internal 
density distributions ranging from r _1 to r~ 3 (e.g., Arquilla & Goldsmith 1985; Chernicaro et al. 1985; Gregorio 
Hctcm et al. 1988; see also Myers 1985). A reasonable mean value for the observationally derived power-law 
distributions is p <~ r~ 2 . The typical sizes of the massive high-density cores are about r c <~ 0.1 pc (see Walmsley 
1995), and the observationally derived core masses are in the range of 10 to 300 M Q (e.g., Sncll et al. 1993). 
The observed properties of ultracompact HII regions (UCHII), on the other hand, indicate that the exciting 
stars (one or several massive stars) are embedded in dense and warm cores, with densities between 10 4 -10 7 cm~ 3 
and temperatures of about 10 2 K (e.g., Churchwell 1990; Cesaroni et al. 1994; Kurtz et al. 1994; Hofncr et al. 
1996; Hurt et al. 1996). The early stages of HII region evolution, then, occur inside these dense cloud cores. 

2.1. HII region expansion at constant densities 

Beginning with Stromgrcn (1939) and Kahn (1954), the expansion and evolution of HII regions has been 
studied with analytical and numerical models (see Yorke 1986; Osterbrock 1989; Franco et al. 1989,1990). For 
a constant photon flux and uniform ambient densities, the evolution has well defined formation and expansion 
phases. During the formation phase, the UV photon field creates an ionization front that moves through the 
gas. Its speed is reduced (by geometrical dilution and recombinations), approaching a value of about twice the 
speed of sound in the ionized gas in, approximately, one recombination time. At this moment, which marks the 
end of the formation phase, the HII region reaches the initial Stromgren size, and the pressure gradient across 
the ionization front begins to drive the expansion of the ionized gas. The expansion is supersonic with respect 
to the surrounding gas and creates a shock wave that accelerates and compresses the ambient shocked medium. 
The ionization front sits behind the shock front during the rest of the evolution, and most of the shocked gas 
is accumulated in the interphase between the two fronts. If either the ionization or the shock front encounters 
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a strong negative density gradient (say, the edge of the cloud) and overruns it, then the HII region enters into 
its champagne phase. Numerical models indicate that there also appears a strong instability that can fragment 
the shocked shell, and may explain the existence of cometary globules and the highly irregular morphologies of 
ionized nebulae (see Garcfa-Segura & Franco 1996). 

2.1.1. The Formation Phase 

Assuming a self-gravitating spherical cloud with a molecular density distribution that includes a central 
core, with radius r c and constant density n c , and an isothermal envelope with a power-law density stratification 
n H 2 { r ) = n c{r/r c ) 2 for r > r c . The total pressure at the core center is 

P(0) = P = ^P 2 A + P(r c ) - |p(r c ) ~ 2 x 1CT 7 ngrg.! dyn cm" 2 , (1) 

where G is the gravitational constant, P(r c ) is the pressure at the core boundary r = r c , n e = n c /10 6 cm -3 , 
and ro.i = r c /0.1 pc. Note that the corresponding molecular mass is 

M c ^^y /2 r 2 ^10 2 P 7 1/2 r 2 , M , (2) 

where P7 = P /10 -7 dyn cm~ 2 . A star located at the cloud center and producing P* ionizing photons per unit 
time creates an spherical HII region. The gas in the HII region is fully ionized and the ion density is simply 
twice the molecular density, m — 2uh 2 ■ The initial Stromgren radius generated by such a star is 



R, 



3 P* 



47r(2n c ) 2 a B 



1/3 

~ 2 x l(T 3 F 4 W 2/3 "o 1/3 pc, (3) 



where as is the hydrogen recombination coefficient to all levels above the ground level, ao — as/2. 6 x 
lCT 13 cm 3 s -1 , P 48 = P»/10 48 s _1 , and n 6 = n c /10 6 cm -3 . Thus, R s < r c and the initial HII region is 
well contained within the core, and the formation phase follows the well known constant density evolution 
described above. Note that for a dusty cloud, the initial radius is even smaller and is given by the 
transcendental equation Rs,dust — Rs e~ T / 3 (Franco et al. 1990), where the optical depth due to dust 

absorption is r = J^ s ' dust cr dust n <ir, and <7d U st is the average dust absorption cross-section per gas particle. 
This approximation agrees to better than 10% with detailed radiative transfer calculations (Diaz et al. 1996). 

The formation phase is completed in a recombination time, and the pressure in the ionized region drives 
a shock into the molecular ambient medium. The HII region now begins its expansion phase. The equilibrium 
temperature in the photoionized region, 7* ~ 10 4 , is achieved in a relatively short time scale. The advance of the 
ionization front is controlled by recombinations in the photoionized gas, and the expansion proceeds in a nearly 
isothermal fashion. Neglecting the external pressure, the main features of the expansion in a constant density 
medium can be derived with the thin shell approximation (sec reviews Ostrikcr & McKee 1988, Bisnovatyi-Kogan 
& Silich 1995, and Garcfa-Segura & Franco 1996). Also, the pressureless expansion in density stratifications can 
be solved with simple approximations to the shock front conditions (Franco et al. 1989,1990). The evolution 
including the external pressure, however, cannot be solved in a closed analytical form, but one can derive limits 
to the main expected features. Obviously, detailed numerical simulations provide an adequate description of 
the evolution in this and more complicated cases. 

2.1.2. The Expansion Phase 

The evolution can be easily derived by assuming the existence of a thin shell with mass M, containing 
all the swept-up ambient gas. This approximation can be applied to HII regions when the fraction of mass 
eroded by the ionization front from the shell is small. This is true for constant, increasing, or mildly decreasing 
density stratifications, but is not applicable for strongly decreasing gradients because the shell is easily eroded by 
photoionization (Franco et al. 1990). Here we consider the constant density case and the thin shell approximation 
can be used without restrictions. Neglecting magnetic fields and self-gravity, the equation of motion of the shell, 
located at a distant R from the central star, is 

4irR 2 (P l -P ) = j t (Mv) 1 (4) 
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Fig. 1. The evolution of the exponent /? as a function of the normalized pressure difference (Pi — P )/Pi- 



where Pi is the internal pressure, Po 1S the external pressure, and v is the shell velocity. Assuming that the 
shell radius can be written as a power-law in time, R = RoC 13 (where Rq is the initial radius of the region, 
C = (*o + t)/t n , and t is a reference initial time), with constant Ro or (3 (i.e., neglecting the existence of terms 
with i? and /3), the right hand side of the equation becomes 



. , /4/3- 



4irp R 2 v 2 . 



(5) 



The assumption of a constant /3, as we will see below, does not hold for cases with P > 0, but the present 
approximations can still be applied in a piece-wise fashion (i.e., one can use them in segments, changing the 
values for Ro and to). 

The equation of motion can then be written as 



Pi - Po = PaV 



(4/3 - 1) 
3/3 : 



and the shell evolution is given by 



dR 
~dt 



3/3 {Pi - P ) 



1 1/2 



.(4/3-1) po 

Thus, one can write the formal solution to the equation of motion simply as 



R — i?o — 



3/3 



(4/3 - l)po 



(Pi - Po) 



1/2 



to dC- 



(6) 



(7) 



(8) 



This formal solution can be applied to HII regions, wind-driven bubbles, and SN remnants (see below). 
Obviously, the integration is not straightforward unless P = 0, or the pressure difference, Pi — P , can be 
written as an explicit function of cither R or t. In general this is not possible, but one can always check the 
behavior at early and late times, when the external pressure can be neglected and when the internal and external 
pressures become comparable. Here we illustrate the behavior in both cases. 



2.1.3. Pressure Equilibrium 



For the case of expanding HII regions, using R$ = R$ as the initial Stromgren radius, we can set the 
solution for the expansion simply as 

R HII (t) = R S C,P . (9) 
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The average ion density inside the HII region at any given time is given by 



< n 4 >= 

and the corresponding internal pressure is 

P 

The region, then, evolves as 
Rhii '- 



3P* 



47T O.B 



3 k 2 T?F, 
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-3/2 
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1/2 
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(4/3 - l) Po 



1/2 



-Pn 



1/2 



(10) 



(11) 



(12) 



At early times, when Pa/Pi << 1, the integration with constant (3 gives Rhii oc C 1 -^/ 4 . Thus, our initial 
assumption for the power-law implies that (3=1 — 3/3/4, and one gets [3 = 4/7. The solution, then, is now 
written as 



Rhii — Rs 



3 k T t 



1/2 



(13) 



Defining a ~ (8 k Ti/S^H) 1 ^ 2 , and recalling the initial definition of the power- law C 4/7 = Rhii/Rs, the 
reference time is simply given by t = (4/7) (Ps jci). Thus, as expected, one recovers the well known law for 
HII region expansion in a pressureless medium with a constant density 



Rhii — Rs ( 1 



7 CiA 



4/7 



Using this explicit time dependence, the internal pressure decreases as 



Pi = P,o 1 + 



7 Cjt 
4i?s 



-6/7 



(14) 



(15) 



where P^o is the pressure at t = 0. 

The expansion continues until Pi — > Po, and the internal pressure tends to a constant value. In this 
limit, the time dependence in the formal solution vanishes, giving 3/3/2 — ► 0. These limits show that, for a 
constant density medium, (3 evolves as a function of the pressure difference from 4/7 to zero. Figure 1 shows the 
evolution of (3, as a function of the normalized pressure difference (Pi — Po)/Pj, for HII regions and wind-driven 
bubbles (Garcfa-Scgura & Franco 1996). The exponents were derived from high-resolution numerical simulations 
performed in one dimension for the evolution inside high density cores. Figure 2 shows the evolution of an HII 
region in a high-density core with P7 = 1. 

When pressure equilibrium is reached, the ion density is simply given by 



^) * 3.6 x 10 4 P 7 T m \ 4 cm- 



(16) 



where P7 = Po/10 7 dyn cm 2 , and Thii,4 = Ti/10 4 K. The equilibrium radius of the H 11 region, then, 
corresponds to a Stromgren radius at this equilibrium density 



Rs 



eq 



2.9 x 10- 2 P 4 V 3 T^ A P~ 



2/3 



pc, 



(17) 



48 o-l 



where P 48 = P*/10 4S s 

For high-pressure cores with r c ~ 0.1 pc, the photoionized regions can reach pressure equilibrium without 
breaking out of the core. The resulting sizes are similar to those of the ultracompact class (see Dyson et al. in 



0.06 



0.04 



R(pc) _ 



0.02 








J I I L 



• Ionization Front 



_L 



_L 



_L 



_±_ 



1 2 ((l .4 .3 4 5 

t (10 yr) 

Fig. 2. The evolution of an HII region, with F48 = 1, in a high-pressure core, with P7 = 1. 



this volume), and indicate that UCHII can be explained by simple pressure equilibrium. This is in agreement 
with the recent results reported by Xie et al. (1996), that show the smaller UCHII are embedded in the higher 
pressure cores. Also note that the equilibrium values with P7 = 1, Thii,4 ~ 1, and F^s ~ 1, are very similar to 
the average sizes and electron densities in UCHII (see Figure 151 of Kurtz et al. 1994). The apparent longevity 
problem of UCHIIs is rooted in the notion that young HII regions should grow fast and reach the expanded 
state on a relatively short time-scale. This statement is false, however, if the external pressure is large and halts 
the expansion at a small radius: in pressure equilibrium, UCHIIs are stable and long lived. 

2. 2. HII Evolution in Decreasing Density Gradients 

If a star is born at a distance smaller than Rs,eq from the core boundary, the speeds of the ionization 
and shock fronts are modified by the negative density gradient. The gas of the HII regions is accelerated in 
supersonic flows and no static solution, in pressure equilibrium exists (e.g., Tenorio-Tagle 1982; Franco et al. 
1989, 1990). Thus, under these conditions the ultracompact stage is indeed a transient phase. Here we assume 
isothermal clouds with p cx r~ 2 , but one can easily find the solutions for the general power-law case, p oc r~ w 
(see Franco et al. 1990). For R s > r c , the initial ionization front reaches the core radius with a speed 



U c ~ 90 a n 3 r 17 



km s 1 , 



(18) 



and in a time scale 



t c ~ 130 a n 3 1 In 



1 - {rJRsY 



(19) 



where m — r c /10 cm. Afterwards, the ionization front enters the density gradient and its speed becomes 



Ui 



Uc 



(R,/r c y - 1 



u(w), 



(20) 



with 



u(w) 



(r c /rif~ w \(R s /r c f + 2wf3 - 3/3 {n/r c ) 1/P } for w ± 3/2 



(r c / n ) 1/2 (R s /r c y- 1-3 In ( n /r c ) 



for w = 3/2, 



(21) 
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whcre ri is the location of the front, and /? = (3 — 2u>) 1 . This defines a critical exponent, corresponding to the 
maximum density gradient that is able to stop the ionization front, 



W f= 2 



Rs 



(22) 



and above which no initial HII radius exists. Note that for R s /r c > 2 the critical value becomes Wf ~ 3/2. 
Obviously, the concept of a critical exponent is not restricted to power-law stratifications and can also be applied 
to other types of density distributions, for instance, exponential, gaussian, and sech 2 profiles (see Franco et al. 
1989). 

For w ^ 3/2, the initial HII region radius can be written as 



R w = g{w)R s , 



with 



- 2w 2w 
~3 + IT 



3-1/3 



2w(3 



(23) 
(24) 



where R s is the Stromgren radius for the density n c . The solution for w — 3/2 is 



R3/2 = r c exp • 



- 1 



(25) 



After the formation phase has been completed in clouds with w < Wf, the HII region begins its expansion 
phase. For simplicity, we assume that the shock evolution starts at t = when R w is achieved. Given that 
the expansion is subsonic with respect to the ionized gas, the density structure inside the HII region can be 
regarded as uniform and its average ion density at time t is 



Pi (t) c ^ \ 2n c )Rl /2 R- m (t), 



(26) 



where [ii is the mass per ion, and R(t) is the radius of the HII region at the time t. For w < 3/2, the radius 
can be approximated by 



R(t) ~ R v 



1 



7-2w 



12 



9 - 4w 



1/2 



Rw 



4/(7-2id) 



(27) 



where Cj is the sound speed in the ionized gas. The ratio of total mass (neutral plus ionized), M s (t), to ionized 
mass, Mi(t), contained within the expanded radius evolves as 



M s (t) 
Mi(t) 



R(t) 

Rid 



(3-2w)/2 



(28) 



This equation indicates: i) for w < 3/2, the interphase between the ionization front and the leading shock 
accumulates neutral gas and its mass grows with time to exceed even the mass of ionized gas, and ii) for 
w = 3/2— w cr it, the two fronts move together without allowing the formation and growth of a neutral interphase. 
Note that the decreasing ratio predicted by equation (28) for w > 3/2 is physically meaningless and it only 
indicates that the ionization front overtakes the shock front (and proceeds to ionize the whole cloud). Thus, 
regardless of the value of the critical exponent for the formation phase, Wf, the expansion phase is characterized 
by a critical exponent with a well defined value, w crit = 3/2, which is independent of the initial conditions. 
Furthermore, this critical exponent w cr i t — 3/2 is not affected by dust absorption (see Franco et al. 1990). 

For 3/2 < w < wj, the ionization front overtakes the shock and the whole cloud becomes ionized. In 
this case, the pressure gradient simply follows the density gradient. The ionized cloud is set into motion, but 
the expanded core (now with a radius identical to the position of the overtaken shock) is the densest region 
and feels the strongest outwards acceleration. Then, superimposed on the general gas expansion there is a 
wave driven by the fast growing core (the wave location defines the size of the expanded core), and the cloud 
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experiences the so-called "champagne" phase. This core expansion tends to accelerate with time and two 
different regimes, separated by w — 3, are apparent: a slow regime with almost constant expansion velocities, 
and a fast regime with strongly accelerating shocks. The slow regime corresponds to 3/2 < w < 3 and the core 
grows approximately as 



r(t) ~r c + 



3 — w 



1/2 



(29) 



where for simplicity the initial radius of the denser part of the cloud has been set equal to r c , the initial size of 
the core. For w — 3 the isothermal growth is approximated by 



r(t) ~ 3.2r c 



Cit 



1.1 



(30) 



For w > 3, the fast regime, the shock acceleration increases with increasing values of the exponent and the core 
expansion is approximated by 



r(t) ~ r c 
where 6 ~ 0.55(w - 3) + 2.8. 



1/2 /S + 2-w\ at 

■'<) V 2 J 77 



2/(5+2— iu) 



(31) 



2.3. Wind-driven Bubbles 



The evolution of the cavity created by a stellar wind, a wind-driven bubble, can also be derived with the 
thin shell approximation described above. The thermalization of the wind creates a hot shocked region enclosed 
by two shocks: a reverse shock that stops the supersonic wind, and an outer shock that penetrates the ambient 
gas. The gas processed by each shock is separated by a contact surface, the contact discontinuity. The kinetic 
energy of the wind is transformed into thermal energy at the reverse shock producing a hot gas {e.g., Weaver et 
al. 1977). For the case of a strong wind evolving in a high-density and dusty molecular core, the properties of 
the cooling are poorly known, but the shocked ambient gas cools down very quickly and a thin external shell is 
formed on time-scales of the order of years. Thus, the ambient gas is collected in a thin shell by the outer shock 
during most of the evolution. The case of the shocked stellar wind is less clear because the cooling time there 
can be substantially longer than in the shocked ambient gas, and it is difficult to define when the thin shell is 
formed behind the reverse shock. Thus, one can simply derive the limits for the evolution of the reverse shock 
in both the adiabatic and radiative modes. 



The density in a steady wind, with a constant mass loss, decreases as 



M 



Pw = 



4irr 2 v r 



(32) 



where M is the stellar mass-loss rate, and is the wind speed. The pressure in the shocked wind region is 
defined by the wind ram pressure, Pwv"^, at the location of the reverse shock and is given by 



MVoo_ 



(33) 



where R rs is the radius of the reverse shock. 



2.3.1. Adiabatic Case 



For an adiabatic bubble evolving in a constant density medium and powered by a constant mechanical 
luminosity, _L W , the thermal energy of the shocked wind region grows linearly with time, E t h = 5 L w t/ 11 
(Weaver et al. 1977). The shocked ambient medium is concentrated in the external thin shell and the the 
bubble radius, R^ is the radius of the contact discontinuity. The thermal pressure of the bubble interior changes 
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as Pi = (5L w t)/(22irRl). This pressure is equal to that given by equation (|3 
and outer shocks are related by mass conservation 



R rs — R 



3/2 



10L W t 



1/2 



and the locations of the reverse 



(34) 



The solution for Rb, then, also provides the evolution of the reverse shock. The initial radius in this case is very 
small, and we simply set Rb = Ro(t/to)@ , where Rq now represents the bubble radius at some reference time to. 
Again, as in the HII region case, the formal solution is 



Rb 



-30 t 30 



(4/3 - l)p 



22ttP 3 



Po 



1/2 



dt. 



(35) 



At early times, Pi >> Po and Po can be neglected in the above equation. The exponent is then defined 
by (3 = 1 + (1 — 3/3)/2, giving the well known adiabatic expansion in a medium with constant density 

Rb oc Ll/ 5 p Q 1 ^ 5 t 3 / 5 . The position of the reverse shock (equation |34|), then, evolves as R rs cx t 2 / b . With this 
time dependence, the internal pressure drops as Pi oc i~ 4 / 5 . At later times, when Pi — > Po, the bubble reaches 
quasi-equilibrium with the ambient gas and the growth is i?b iGq cx (L w / Po) 1 ' 3 * 1 ' 3 . In pressure equilibrium, the 
radius of an adiabatic bubble grows at a slow rate, but no steady state solution exists during this stage. The 
final radius of the reverse shock is simply given by the balance between the external and the wind ram pressures. 



The wind density at equilibrium is p v 



Pq I v 2 ^ , and the location of the reverse shock is 



R 



rs,cq 



Mvo_ 
4ttP 



"I 1/2 



2.3 x 1CT 2 



M 6 Woo, 8 



1/2 



pc, 



(36) 



where Me = M/10 6 M© yr x , and 1*00,8 = Voo/10 8 cm s 1 . Using mass conservation in the shocked wind 
region, M sw — Mi, the quasi-equilibrium radius of an ultracompact wind-driven bubble is given by 



Rb.cq — Rrs,eq I 1 



1/3 



87rP P 3 S!Cq/ 

Clearly, for adiabatic bubbles, j3 goes from 3/5 at early times to 1/3 at late times. 

2.3.2. Radiative Case 



(37) 



Once radiative losses become important, the hot gas looses its pressure and the bubble collapses into a 
simple structure: the free-expanding wind collides with the cold shell and the gas is thermalized and cools down 
to low temperatures in a cooling length. At this moment, the shell becomes static at the radius P rS)eq - 

If the bubble becomes radiative before the final radius, R rSjeq , is reached, the shell is pushed directly by 
the wind pressure. For this radiative bubble case, the formal solution is 



Rb = 



33 



MVoct'o 



2ii 



(4/3- 



l)po I AkRIW 



Pn 



1/2 



dt. 



(38) 



The early times solution, with Pi » Pq, gives the relation /3 = 1 — /3, and one recovers the well known solution 
for radiative bubbles Rb cx t 1 / 2 {e.g., Steigman et al. 1975). The ram pressure now evolves as Pi cx t~ l . At late 
times, one gets /3 — > and the shell reaches the final equilibrium radius (eqn. |36|). Thus, the exponent in this 
case evolves from 1/2 to 0. Note that if the bubble starts on an adiabatic track and becomes radiative before 
pressure equilibrium is achieved, the exponent varies from 3/5 to 1/2, and then to zero. Figure 3 illustrates the 
evolution of an ultracompact wind-driven shell in a high-pressure core. 
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Fig. 3. The evolution of an ultracompact wind-driven bubble, with Mq — 1 and v Q 
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3. DISCUSSION 
3.1. Young Stars 

The history of the pressure in a star forming region, then, follows a somehow simple scheme. The initial 
properties and pressure of the gas in a star forming cloud is defined by self-gravity. Once young stars appear, 
their energy input modifies the structure and evolution of the cloud. This is particularly true for massive stars, 
their radiative and mechanical energy inputs are even able to disrupt their parental clouds. In the case of 
the dense cloud cores, the sizes of either HII regions or wind-driven bubbles are severely reduced by the large 
ambient pressure (Garcia-Scgura & Franco 1996). In fact, the pressure equilibrium radii of ultra-compact HII 
regions are actually indistinguishable from those of ultra-compact wind-driven bubbles, and they could be stable 
and long lived. Actually, Xie et al. (1996) have recently found evidence indicating that the smaller UCHII 
seem to be embedded in the higher pressure cores. 

The situation is completely different when the stars are located near the edge of the cloud core. The 
resulting HII regions (and also wind-driven bubbles) generate supersonic outflows. Cases with w > w cr i t lead 
to the champagne phase: once the cloud is fully ionized, the expansion becomes supersonic. For spherical 
clouds with a small constant-density core and a power-law density distribution, r~ w , outside the core, there is a 
critical exponent (w cr it = 3/2) above which the cloud becomes completely ionized. This represents an efficient 
mechanism for cloud destruction and, once the parental molecular cloud is completely ionized, can limit the 
number of massive stars and the star formation rate (Franco et al. 1994). For a cloud of mass MqmCi with 
only 10% of this mass concentrated in star-forming dense cores, the number of newly formed OB stars required 
for complete cloud destruction is 



AT or, MgmC&H ,„ q , 

^48 (c 4 ,15^MS,7) 6/5 

where M GMC 5 = M c /10 5 M Q , n 3 = n /10 3 cm~ 3 , c iil5 = Cj/15 km s _1 , and tMS,7 is the main sequence lifetime 
in units of 10^ yr. This corresponds to a total star forming efficiency of about ~ 5 % (larger average densities 
and cloud masses can result in higher star formation efficiencies). 

Summarizing, photoionization from OB stars can destroy the parental cloud in relatively short time scales, 
and defines the limiting number of newly formed stars. The fastest and most effective destruction mechanism is 
due to peripheral, blister, HII regions, and they can limit the star forming efficiency at galactic scales. Internal 
HII regions at high cloud pressures, on the other hand, result in large star forming efficiencies and they may be 
the main limiting mechanism in star forming bursts and at early galactic evolutionary stages (see Cox 1983). 
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Fig. 4. Wind-driven bubble from an evolved 35 M star in a high pressure medium, a) Left scale: terminal 
wind velocity (km s _1 ); right scale: mass loss rate (M Q yr -1 ); horizontal scale: time (millions of yr). (b)-(d): 
Evolution of the wind-driven bubble. The gas density (cm -3 ) is plotted in a logarithmic scale versus the radial 
distance (pc). Evolutionary times (upper- left corner of each panel) are given in million years. 

3. 2. Slow Winds from Evolved Stars 

The evolution of evolved stellar winds in high-pressure regions is discussed in Franco et al. (1996), and here 
we just repeat the relevant parts. As the cloud is dispersed, the average gas density decreases and the newly 
formed cluster becomes visible. The individual HII regions merge into a single photo-ionized structure and the 
whole cluster now powers an extended, low density, HII region. The stellar wind bubbles now can grow to larger 
sizes and some of them begin to interact. As more winds collide, the region gets pressurized by interacting 
winds and the general structure of the gas in the cluster is now defined by this mass and energy input (Franco 
et al. 1996). 

Given a total number of massive stars in the cluster, Nob, and their average mass input rate, < M >, the 
pressure due to interacting adiabatic winds is 

N OB <M> Cl in _ 8 Ar 2 < M 6 > c 20 oo , _ 2 , An ^ 
Pi ~ — , ~ 10 9 d y n cm > ( 40 ) 

clus pc 

where r pc — r c i us /l pc is the stellar group radius, N 2 = Aos/lO 2 , < M 6 >=< M > /10 -6 M yr -1 , and 
C2000 = Ci/2000 km s -1 is the sound speed in the interacting wind region. This is the central pressure driving 
the expansion of the resulting superbubble before the supernova explosion stage. For modest stellar groups 
with relatively extended sizes, like most OB associations in our Galaxy, the resulting pressure is only slightly 
above the ISM pressure (i.e., for N 2 ~ 0.5 and r pc ~ 20, the value is Pi ~ 10~ n dyn cm -2 ). For the case 
of rich and compact groups, as those generated in a starburst, the pressures can reach very large values. For 
instance, for the approximate cluster properties in starbursts described by Ho (1996), r pc ~ 3 and N 2 > 10, the 
resulting pressures can reach values of the order of Pi ~ 10~ 7 dyn cm~ 2 , similar to those due to self-gravity in 
star forming cores. At these high pressures, the wind of a red giant (or supergiant) cannot expand much and 
the bubble reaches pressure equilibrium at a relatively small distance from the evolving star. Thus, the large 
mass lost during the slow red giant wind phase is concentrated in a dense circumstellar shell. Figure 4 shows 
the evolution of a wind-driven bubble around a 35 M Q star. Fig. 4a shows the wind velocity and mass-loss 
rate (dashed and solid lines, respectively: Garcfa-Segura, Langer & Mac Low 1996). The simulations are done 
only over the time spanning the red supergiant and Wolf-Rayet phases, and assume that the region is already 
pressurized by the main sequence winds from massive stars. 

We have used the AMRA code, as described by Plewa & Rozyczka (1996). During the RSG phase the 
wind-driven shell is located very close (R w 0.04 pc) to the star due to a very low wind ram-pressure from the 
wind (Fig. 4b). Later on (Fig. 4c), the powerful WR wind pushes the shell away from the star to the maximum 
distance of R w 0.54 pc. Still later, when the wind has variations, the shell adjusts its position accordingly, 
and reaches the distance R w 0.3 pc at the end of simulation (Fig. 4d). It must be stressed that the series 
of successive accelerations and decelerations of the shell motion during the WR phase will certainly drive flow 
instabilities and cause deviations from the sphericity assumed in our model. The role of these multidimensional 
instabilities in the evolution of the shell is currently under study (with 2-D and 3-D models), and the results will 
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be presented in a future communication. Regardless of the possible shell fragmentation, however, when the star 
explodes as a supernova, the ejecta collides with a dense circumstellar shell. This interaction generates a bright 
and compact supernova remnant, with a powerful photoionizing emission (i.e., Terlevich et al. 1992; Franco et 
al. 1993; Plewa & Rozyczka 1996), that may also be a very strong radio source, like SN 1993J (see Marcaide 
et al. 1995). If the shell is fragmented, the ejecta- fragment interactions will occur during a series of different 
time intervals, leading to a natural variability in the emission at almost any wavelength (see Cid-Fernandes et 
al. 1996). This type of interaction is also currently under investigation, and further modeling will shed more 
light on the evolution of SN remnants in high-pressure environs. 

JF and GGS acknowledge partial support from DGAPA-UNAM grant IN105894, CONACyT grants 400354- 
5-4843E and 400354-5-0639PE, and a R&D Cray research grant. The work of TP was partially supported by the 
grant KBN 2P-304-017-07 from the Polish Committee for Scientific Research. The simulations were performed 
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Institut fiir Astrophysik. 
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